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ABSTRACT 

When modeling the three-dimensional hydrodynamics of interstellar material 
rotating in a galactic gravitational potential, it is useful to have an analytic ex- 
pression for gravitational perturbations due to stellar spiral arms. We present 
such an expression for which changes in the assumed characteristics of the arms 
can be made easily and the sensitivity of the hydrodynamics to those character- 
istics examined. This analytic expression also makes it easy to rotate the force 
field at the pattern angular velocity with little overhead on the calculations. 

Subject headings: hydrodynamics — Galaxy:structure — galaxies:spiral, struc- 
ture 



1. INTRODUCTION 

In this paper we present analytic expressions for the perturbation of the galactic ax- 
isymmetric gravitational potential due to redistribution of part of the stellar disk mass into 
spiral arms, and for the distribution of density responsible for that potential. Adjustable 
parameters include the number of arms, N, the pitch angle, a, the radial scale length of the 
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dropoff in density amplitude of the arms, i? s , the midplane arm density, p Q at fiducial radius 
r Q , and the scale height of the stellar arm perturbation, H . 

The amplitude of the spiral density distribution whose gravitational potential we set 
out to find is given by: 

p A (r, z) = Po exp ^—j^J sech 2 {^j ■ (1) 

Modulating this by a simple sinusoidal pattern in 0, following a logarithmic spiral with a 
pitch angle a, the overall density perturbation is 

p(r, 0, z) = p A (r, z) cos(7) (2) 

where 



7 = iV 



\n(r/r ) 

- M r o) - 



tan(a) 



(3) 



More complicated azimuthal arm structures can be constructed with linear combinations 
of these solutions of the form 

P(r, 0, z) = p A (r, z) C n cos(n7). (4) 

n 

In a particularly interesting example, the density behaves approximately as a cosine squared 
in the arms but is separated by a flat interarm region occupying half the volume. It has 
three terms in its sum, with C\ = 8/(37r), C2 = 1/2, and C3 = 8/ (157r). The resulting phase 
pattern is compared with that of a simple sinusoid in Figure 1. 

An important feature of such a perturbation is that its average density is zero. Its po- 
tential can thus be added to observationally constrained models for the azimuthally averaged 
potential without altering the latter. In addition, because the assumed arm perturbation 
extends over all radii, it is important that the average density be zero at both large and small 
radii where the arms do not actually exist. The radial exponential damping introduced in 
Equation 1 is also useful in this regard. Thus the gradient of the perturbation potential in 
the calculation region is provided predominantly by the local distribution of material. 



2. THE POTENTIAL 

In developing an expression for the potential, we first considered an infinite sinusoidally 
oscillating density pattern in a rectilinear coordinate system. Mathematically it was equiv- 
alent to Equations 1 and 2 with R s infinite, 7 replaced with kx, and y extending from + 
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to - oo, normal to the wave vector. For this density distribution, we evaluated the potential 
numerically and fitted the results versus phase (kx) and z with a simple functional form. 
We then remapped the solution's phase into the desired spiral pattern and reintroduced the 
radial dropoff factor, exp(— r ^ r2 -), obtaining a trial potential function. 



The resulting potential has three functional parameters dependent on radius: 

K„ . " N 



r sin (a) 
f3 n = K n H(l + 0AK n H) 

l + K n H + 0.3(K n H) 2 



D 



1 + 0.3fT n iJ 



(5) 
(6) 
(7) 



With these parameters, the trial potential corresponding approximately to the density dis- 
tribution of Equation 4 is 



$(r, 0, z) = -4%GHp exp 



r — r c 



E 

n 



C n 
K n D n 



cos(n7) 



sech 



1 /9n 



(8) 



THE DENSITY 



Our next step was just an inversion of attitude, from considering the density function 
of equations 1 through 4 as primary and the potential function of equation 8 as an approx- 
imate solution, to the inverse perspective. Equation 8 is our final equation for the potential 
perturbation, and by evaluating 

we find the exact corresponding density function. Because of the radial dependence of the 
parameters K n , /3 n , and D Q , the radial derivatives of the Laplacian are particularly messy. 
The full solution is given in the Appendix; for each n there is one dominant term for small 
H/r and not too large z/H such that: 



p(r, 0, z) w p exp 



r -r c 



Pn 



cos(n7) 



sech 



K n z 



(2+AO 



(10) 



Of course, for particular choices of parameters, the Laplacian can be evaluated numerically to 
determine the corresponding density function. We did so to check the exact density solution 
of the Appendix, for example. 
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4. EXAMPLES 

We next present the results for two cases, one with simple sinusoidal arms and one with 
the more concentrated arms, as per Figure 1. The parameters in both cases are N = 2, 
a = 15°, R s = 7 kpc, p Q = mn , n Q = 1 atom cm' 3 at r G = 8 kpc, and H = 0.18 kpc (chosen 
to match the scaleheight of the thin stellar disk of Dehnen and Binney, 1998, Model 2). The 
average mass per atom is m = (14/1 l)m#. Unless otherwise noted, results are shown versus 
r and z along a radial cut passing through maximum arm density at r = 8 kpc. 

Figures 2 and 3 show the perturbing gravitational potential and corresponding density 
found from the Laplacian of the potential, for the sinusoidal arm case. Figures 4 and 5 
are similar but for the concentrated arms case. In both cases, the density functions are 
almost indistinguishable from the assumed function of Equation 4 or the dominant term 
form of Equation 10; the differences are much smaller than the uncertainty in how the true 
arms should be represented. As a result, Equation 4 (or 10) can be used with considerable 
confidence as an approximation to the density distribution responsible for the potential. 

It is somewhat easier to compare these two cases by examining separately their radial and 
vertical behaviors. Figures 6 and 7 compare the midplane densities and potentials versus 
radius. The maxima and minima of the density functions are bounded by the decaying 
exponential envelope with R s = 7 kpc. This value for the radial dropoff scale was chosen 
so that the amplitude of the potential variation would not depend strongly on radius. Our 
"sinusoidal" case is sinusoidal with phase (and in 0), but in the logarithmic spiral, the 
effective wavelength is proportional to radius. The "concentrated" case has the desired flatter 
interarm density and sharper peaks at the arms. The difference is moderated somewhat in 
the potential function, but the flatter peaks (interarm) and sharper valleys of the arms are 
apparent in the concentrated arm case. Figures 8 and 9 compare the vertical structures at 
the location of an arm. The greater density and deeper and sharper potential well of the 
concentrated case are evident. 

Comparison of Figures 8 and 9 emphasizes the much greater scale height of the po- 
tentials, compared to that of the responsible densities, evident also in Figures 2 through 
5. The potential scale height depends on both the density scale height and on the radial 
wavelength. For an infinite disk (infinite radial wavelength), the vertical scale height is also 
infinite, having a constant gradient outside the mass distribution. Our results successfully 
model the useful cases with K n H small to moderate. 

The impressions invoked by the density distributions of Figures 3 and 5 can be somewhat 
misleading. These densities must be regarded as perturbations to an azimuthally uniform 
stellar disk with the same vertical scale height. In Figure 10, a disk component with the 
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same radial dropoff and scale height as the perturbation, and just sufficient amplitude to 
make the net density everywhere positive, has been added to the perturbation density. In 
the figure captions, this disk is referred to as the Net Mass Disk, as it provides the material 
for the arms. This picture provides a better idea of the assumed density structure of the 
arms. 

Those tempted to find the potential of the above combination of disk plus perturbation 
may wish to consider an alternative approach to evaluation of the spiral arm potential that 
has been pursued by Martos, Pichardo, Moreno, and Espresate (2002); they specify a positive 
density pattern and evaluate the potential with mixed analytic and numerical techniques. 

Having created this intuitive view of the density distribution, it is used in Figure 11 to 
illustrate the qualitative distribution of arm material in the galactic midplane. 

In Figures 12 and 13, various amounts of perturbation density are shown added to a 
representative full stellar disk. The disk parameters were chosen to approximate the thin 
stellar disk in Model 2 of Dehnen and Binney (1998). It has a sech 2 (z / H) form with a scale 
height of 0.18 kpc, a radial dropoff scale of 2.4 kpc, and an atomic number density of 3.21 
cm~ 3 at r = 8 kpc. In the sinusoidal arm case, the added perturbation has, at 8 kpc, an 
amplitude equal to the fractions of the total 3.21 cm -3 disk density that are shown in the 
legend. A fraction of 0.4, for example implies total interarm and arm densities of 0.6 and 1.4 
times the unperturbed density. At smaller radii, the perturbation is a smaller fraction of the 
full disk because of their different radial dropoff scales. 1 In the concentrated arm case, the 
normalization is slightly different to provide similar arm to interarm density difference at 
fixed radius, as per Figure 1. In this case, a 0.4 perturbation fraction implies total interarm 
and arm densities of 0.8 and 1.6 times the unperturbed density at r = 8 kpc. (The average 
density is the same in the two cases because the concentrated arms are narrower.) These 
figures provide a sense of the degree of modulation required for an observable density contrast 
in the presence of the steep gradient of the unperturbed disk. Rather large perturbations 
are required just to flatten the density profile. When the perturbation is large enough to 
give the total density a local maximum, that maximum is shifted radially inward from the 
perturbation peak. 

Figures 14 and 15 further examine the stellar density distribution. Figure 14 shows the 



1 The extreme difference between the 2.4 kpc radiaf dropoff scafe for the average thin disk density and 
our choice of 7 kpc for the radial dropoff scale of the perturbation density was not guided by observations. 
We are using the perturbation in 3D MHD modeling of the response of the interstellar medium to the stellar 
arm potential, and did not want to introduce, at least initially, a strong radial gradient in the depth of the 
arm potential well. 



- 6- 



vertical distribution of the density versus radius for the stellar disk described above plus 
"0.57" fraction sinusoidal arms. Figure 15 is similar, but the perturbation is shown added 
to the entire axisymmetric galactic density distribution of Model 2 of Dehnen and Binney 
(1998). (Actually, the potentials were added and then the density evaluated by taking the 
Laplacian numerically.) In the spiral arm modeling of Gomez and Cox (2002), these are the 
parameters used for the two-armed case. It had been our prejudice that a 57% modulation of 
the disk density at r = 8 kpc represented a rather extreme situation, and we were therefore 
surprised to learn its rather modest effect on the density structure in comparison to the 
radial gradient. 

These results can be compared with measurements of the actual red light profiles of 
spiral galaxies, if the contamination due to red supergiants can be neglected, to provide a 
reasonable estimate for the arm perturbation amplitude. This was done, for example, for M51 
by Rix and Rieke (1993) and for several other galaxies by Rix and Zaritsky (1995). They 
concluded that in galaxies with strong arms, the arm to interarm stellar density contrast 
was roughly a factor of two. They further found that in some two-armed spirals, the non- 
axisymmetric structure was dominated by the fundamental sinusoidal component, implying 
very broad arms (as in our "sinusoidal" case); but in cases with very strong arms there was a 
significant higher harmonic, implying narrower arms (as in our "concentrated" case). Thus, 
density contrasts at the upper end of those we have explored seem to be relevant to at least 
strong armed galaxies, as do both varieties of the arm modulation functions of Figure 1. 



ROTATING THE PATTERN 



If one is calculating the hydrodynamics in the inertial frame, it is necessary to rotate the 
perturbing potential at the pattern angular velocity Q. For <fi increasing counterclockwise 
and a clockwise rotation, examination of Equations 3 and 8 reveals that each term in the 
sum consists of an amplitude function 



A n (r, z) = —4iiGHp exp 



K n D n 



sech 



K n z 



multiplied by a phase function cos(nr), where 



T(r,0,t) =-f + NQt = N 



(j) + Vtt - p (r o ) - 



In (r/r D ) 



(11) 



(12) 



tan(aj 

Thus, for each term in the sum of Equation 8, only the function A n (r, z) needs to be evaluated 
and stored. At any time, they can be combined according to 



$(r, 0, z,t) = ^ A n (r, z) cos (nF(r, 0, t)) 



(13) 
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THE EXACT DENSITY DISTRIBUTION 



The exact density distribution reflected in the nth term of the potential of Equation 8 
is given by: 



Pn(r,<j),z) = C n p 



where 



and 
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Pn 
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sech 



Pn 



0A(K n Hf 



+1.2(K n Hf 



Pn 



K n z 

Pn 



sech 



tanh 



K n z 

Pn 
KnZ 

Pn 



'Pn 



We have stated in the text that, except at very small r or large z, the terms other than 
that in Equation 10 are very small. Trying to state this more carefully, we note that for 
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K n H ^> 1, E and rE' approach K n z and —K n z, respectively. In that limit, the density 
function approaches 



the third term will exceed the first. This occurs even with our standard parameters; its effect 
can be clearly seen in Figures 3, 5, and fO, as a shift in phase of the density perturbation 
above about 500 pc. So, it is inaccurate to say that the small terms are everywhere negligible 
compared to the leading term. What is true is that the small terms will dominate at high z 
and small r, but in practical cases of moderate z/r and K n H y£ 1, the small terms will be 
small where the density is large, and not very noticeable where it is small. They drop off 
slowly with z, about like the perturbation potential. 
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Fig. 1. — Comparison of Modulation Functions for Sinusoidal and Concentrated Arms. 
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Fig. 2. — Gravitational Potential for Sinusoidal Arms. This plot and all subsequent ones 
versus r and z are for a radial cut at constant which intersects an arm at 8 kpc. Parameters 
are given in Section 4. 
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Fig. 3. — Distribution of Density Perturbation for Sinusoidal Arms. Corresponding potential 
is shown in Figure 2. 
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Fig. 4. — Gravitational Potential for Concentrated Arms. 
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Fig. 5. — Distribution of Density Perturbation for Concentrated Arms. Corresponding 
potential is shown in Figure 4. 
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Fig. 8. — Densities of the Two Cases at r = 8 kpc, Versus Height. 



-17- 




Fig. 9. — Potentials of the Two Cases at r = 8 kpc, Versus Height. 
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Fig. 10. — Density Distribution of Concentrated Arm Perturbation Plus Net Mass Disk. See 
text for description of this intuitive distribution of the arm density, rather than the arm 
density perturbation, which has average density zero. 




Fig. 11. — Midplane Arm Density for the Concentrated Arm Case. The domain is 22 kpc 
on a side. The rapidly oscillating density in the inner 3 kpc is not shown. Density includes 
net mass disk with the same scale height, radial dropoff rate, and total mass as the arms. 
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Fig. 12. — Midplane Density of Disk+Arms for Various Arm Amplitudes, Sinusoidal Arm 
Case. The legend shows the fractional amplitude relative to the unperturbed disk at a radius 
of 8 kpc. 
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Fig. 13. — Midplane Density Versus Radius of Disk+Arms for Various Arm Amplitudes, 
Concentrated Arm Case. See text for amplitude description. 
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Fig. 14. — Density Distribution of Disk Plus Arms for 0.57 Amplitude, Sinusoidal Arm Case. 
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Fig. 15. — Full Model Density. Figure shows the full density represented by Model 2 of 
Dehnen and Binney (1998) plus a sinusoidal arm perturbation with our standard parameters 
(see Section 4) with an amplitude at 8 kpc of 57% of the thin disk density. This is the two- 
armed spiral model used by Gomez and Cox (2002) to explore the gaseous response in MHD. 
Note that the arm/interarm contrast is small even at this amplitude. 



